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Abstract 

Recently I. Tokatly and O. Pankratov ("TP", Phys. Rev. B 60, 15550 (1999) 
) used velocity moments of a classical kinetic equation to derive a hydrody- 
namic description of electron motion in a degenerate electron gas. Indepen- 
dently, the present authors (Theochem 501-502, 327 (2000)) used considera- 
tions arising from the Harmonic Potential Theorem (Phys. Rev. Lett. 73, 
2244 (1994)) to generate a new form of high-frequency hydrodynamics for in- 
homogeneous degenerate electron gases (HPT-N3 hydrodynamics). We show 
here that TP hydrodynamics yields HPT-N3 hydrodynamics when linearized 
about a Thomas-Fermi groundstate with one-dimensional spatial inhomno- 
geneity. 

67.55.Fa, 72.30.+q, 73.21.-b, 52.30.-q 
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I. INTRODUCTION 



The prediction of collective oscillations in degenerate electron gases is a very old problem 
that has been addressed over the years at many different levels of theory. In the limit of 
an infinite homogeneous electron gas it has been possible to apply relatively sophisticated 
microscopic theory beyond the Random Phase (RPA) approximation0'i. The main current 
issue in the degenerate case of the uniform gas is the effect of electron-electron correla- 
tions on the relevant response functions, leading to modified dispersion relations at finite 
wavenumber, possibly presaging quantum phase transitions. Interestingly, however, the ho- 
mogeneous plasmons both in two and three dimensions are completely determined by sum 
rules in the zero- wavenumber limit, so that elementary hydrodynamic arguments give the 
correct answers in this limit. 

The situation is not so well advanced for inhomogeneous degenerate electron gases, where 
high-level dynamic microscopic theory is unwieldy, though the RPA has been carried through 
some time ago for simple-metal surface geometries!. Calculations of the RPA type are now 
possible at the numerical level for molecules and possibly larger systems via a number of 
quantum-chemistry-inspired packages. Some form of Time Dependent Density Functional 
Theory can include the correlation physics approximatelyi'S. All of these microscopic ap- 
proaches are rather computationally intensive, however. 

In order to deal with long-wavelength response in a computationally simple fashion, 
Blochi and others very long ago introduced a hydrodynamic style of treatment for degener- 
ate electron systems. Various workers@'i'i'0'0' have used this idea in linearized form with 
the assumption that the pressure perturbation is proportional to the density perturbation. 
In the early treatments, inhomogeneity of the gas was envisaged only in terms of edges 
bounding an otherwise uniform gas, and these edges were treated via postulated boundary 
conditions. Such approaches, often with a "hard-wall" boundary conditionS, have given a 
largely reasonable description of collective plasmon modes in confined systems such as thin 
metal slabs and quantum wells, with however some shortcomings to be discussed below. Be- 
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cause plasmons are high-frequency phenomena, the connection between pressure and density 
perturbations in a gas of uniform density n, 

5p = mP'^Sn (1) 

had to be used with a nonstandard "high-frequency" value of the coefficient f3 derivable from 
microscopic Lindhard response theory, namely 

(3^ = hUn) = h''{37r'ny/ym^ (2) 





instead of the low-frequency value 



= IvUn) (3) 



derivable from the static pressure-density relation for a degenerate electron gas. 

An obvious defect of the hydrodynamic approach is that it can be guaranteed only to 
describe the "classical" size effects such as standing plasmons in which a finite number 
of nodes of the electron density perturbation fit across the thin dimension of an electron 
gas layer. By contrast, "quantum" size effects involve details of wavefunctions and are 
not necessarily given by hydrodynamic theory: these size effects are related directly to 
transitions between discrete quantum states due to spatial confinement of the electrons. 
These two categories of mode are not fully distinct, howeverEl: for e xample, the sloshing or 
Kohn mode of the electron gas in a parabolic quantum well can be regarded either as due 
to transitions between two Harmonic Oscillator states, or as an odd combination of surface 
plasmonsEl. Nevertheless, it is clear that in general there can be many more modes in a 
microscopic treatment than in the corresponding hydrodynamic one. 

Apart from the lack of the complete spectrum of quantum size-effect modes, the hy- 
drodynamic approach, as usually applied, has other drawbacks for confined electron gases. 
Firstly, it has mostly been used under the assumption that the electron gas is spatially 
homogeneous. This assumption then requires a separate treatment of edge effects, via some- 
what arbitrary boundary conditions at the edge of the uniform-gas region. For example, the 



hard-wall boundary conditions and uniform spatial profile frequently used to analyze plas- 
mons on charge-neutral electron gases of finite width0 give a reasonable description of the 
2D plasmon and sloshing modes (equivalent to even and odd combinations of surface plas- 
mons) but not the multipole surface plasmons, which require a selvage of finite width. For 
an electron gas under parabolic confinement (non-neutral gas) on the other hand, the same 
hard-wall condition does not seem to be so appropriate0Jii, and "free" boundary conditions 
have been suggested along with some further assumptions.. 

The hydrodynamic approach is particularly tricky for plasmon modes that are intrin- 
sically tied to inhomogeneity, such as the multipole surface plasmons of a neutral jellium 
surface. These modes depend intimately on the detailed decay of the electron density at the 
edge, and this cannot be treated by a simple, single boundary condition on a uniform gas. 
There have been attempts within hydrodynamics to model the edge by a series of steps, 
each with its own distinct but internally homogeneous densityEl, in which case additional 
boundary conditions must be used. HydrodjTiamic models have also been investigated with 
a smoother but somewhat arbitrary surface density profilJ'0'@. The problem then is that 
the profile is not internally generated, and in some of these approaches the pressure term 
was not properly consistent with the chosen inhomogeneous density profile. 

In connection with inhomogeneous or bounded degenerate gases, we have taken the point 
of viewlli'0 that hydrodynamics should be applied to a continuously varying model of edges 
(or of inhomogeneity in general). In addition, one should use a smooth solution of the 
groundstate problem that is consistent with the hydrodynamic approximation chosen. Then 
arbitrary boundary conditions can be kept to a minimum or removed altogether. In ad- 
dition, the pressure term should be properly adapted to the inhomogeneity. For plasmon 
applications on a finite slab of electron gas, we tried the obvious inhomogeneous generaliza- 
tion of the high-frequency Bloch approach by using a high-frequency inhomogeneous pres- 
sure coefficient (0) /5^(no(r)) = |f|^(no(r)), where no{f) is the selfconsistent Thomas- Fermi 
groundstate density profile. This approach was foundlil, however, to violate the Harmonic 
Potential Theorem (HPT). This theorem is an extension of the Generalized Kohn Theorem^ 



and states among other things that, for systems confined by a harmonic external potential 
(e.g. the electron gas in a parabolic quantum well), there exists a plasmon mode in which 
the inhomogeneous groundstate density (and groundstate many-body wavefunction) oscil- 
lates rigidly at the bare harmonic oscillator frequency. This violation was subsequently^ 
discussed in a more general fashion, from the viewpoint of invariance under transformation 
to an accelerated reference frame@. In Eqs. (21)- (24) of the original HPT pap a mod- 
ification of coUisionless Bloch hydrodynamics was suggested in order to satisfy the HPT. 
That version of hydrodynamics was later found to violate Newton's third law (momentum 
conservation), but this was rectified to produce a theory that we here term the HPT-N3 
hydrodynamics^. HPT-N3 theory was applied successfuUyil to describe both the Kohn and 
standing-plasmon modes of a parabolic quantum well gas. We have also been able§ to use 
HPT-N3 hydrodynamics to see multipole surface plasmons at the surface of a near-neutral 
electron gas, without using additional boundary conditions. 

It is this modified HPT-N3 hydrodynamic theory that will be discussed further here and 
compared with another hydrodynamic theory. Specifically, Tokatly and Pankratov ("TP") 
in two recent papers ("TPl"ii, "TP2"lli) derived a hydrodynamic theory valid both for 
low and high frequencies, by taking moments of an assumed semiclassical kinetic equation. 
They did not, however, explicitly examine the case of spatially inhomogeneous groundstates 
(other than via assumed boundary conditions on a homogeneous gas). We will show here 
that this TP theory agrees with the HPT-N3 theory in the coUisionless limit, when taken to 
second order in gradients and linearized about the appropriate Thomas-Fermi groundstate 
of an electron gas with general one-dimensional inhomogeneity. Since the two derivations of 
the pressure term are quite different, this agreement tends to support the validity of both 
theories for the case of ID inhomogeneity. 

The present work is organized as follows. In Section II a brief summary is given of the 
HPT-N3 theory. In section HI the coUisionless TP equations through second order are quoted 
for the simplified case of ID spatial variations. In Section IV these ID TP equations are 
linearized about the inhomogeneous TF groundstate and the result is shown to be precisely 



the HPT-N3 theoryE3. Section V is devoted to discussion and conclusions. 



II. SUMMARY OF LINEARIZED HPT-N3 HYDRODYNAMICS 

The exact equation of motion (Euler, Newton II) for the fluid velocity u{f, t) is 

d ^ - - 

mn(— + u.V)u = -V.P + nF (4) 
ot 

where P(r, t) is the stress-pressure tensor (not known exactly in general), n{f, t) is the num- 
ber density, and F (r, t) is the (external plus internal) force per electron due to interactions, 
usually treated via some form of mean-field theory. The pressure force per electron, i.e. the 
force due to streaming of particles from one fluid element into another, is F^^^^^^ = — V.P. 
The form (^) guarantees that the net force on the system due to pressure is zero: specifically, 
the total system pressure force is 

Linearizing (^) about a non-streaming [uq = 0) equilibrium state we obtain 

9 - ^ -* 

Tnno{r)—Ui{f,t) = -V ■ Pi(r, t) + ni(r, t)Fo(r) + no(r)Fi(f, t) (5) 

The approach introduced previously^ in order to accommodate the HPT in a Bloch-style 
inhomogeneous hydrodynamics was to use the fluid displacement x defined by 

x{r,t) = [ u{f,t')dt' (6) 

J to 

to write the density perturbation rii as the sum of a compressive part uia and a displacive 
part tiib: 

ni(f, t) = - V ■ (no (f) f (f, t)) (7) 

= niA{r,t) + niB{r,t) (8) 

niAir, t) = -no(r)V ■ f (f, t) (9) 

niB^r, t) = -f (f, t) ■ Vno(r) (10) 
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Equation (|^) is exact to linear order and can be obtained from a time integration of the 
linearized continuity equation (conservation of electrons). Eqs. are also exact. Note 

that riiB represents the density perturbation that would result if the groundstate density 
distribution were locally displaced rigidly without compression, while Misrepresents the 
density perturbation due to compression. For Kohn-mode motion in which the groundstate 
is translated rigidly, Uib = ni while Hia = because there is no compression. By contrast, 
for oscillations about a spatially homogeneous equilibrium state for which Vno = 0, we have 
niA = ni while tiib = 0. 

Thus in order to satisfy the known Lindhard limit for high-frequency response of a uni- 
form gas, it was proposedEl that riiA be associated in general with the high-frequency pressure 
coefficient /?^(oo) = |t>|,. Similarly to obtain the correct HPT behavior for inhomogeneous 
harmonically-confined systems it was further proposedlll that Uib be associated with the 
low-frequency coefficient /?^(0) = ^Vp , because the static groundstate pressure profile must 
be rigidly translated in Kohn-mode motion. Then the pressure perturbation (assumed0'0 
to be a scalar rather than the more general tensor) is 

pf™(f,t) = mvUnom%A{r,t) + \n,B{r,t)] (11) 
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Pu, = 5,,J>f ^^^^ (12) 

where 

4(no(r)) = ffm'\3-K^no{r)f/'' (13) 

is the square of the local Fermi velocity corresponding to the groundstate number density 
no (r). Although motivated by the limiting uniform and quadratically-confined cases, Eq. 
(|n|) was proposed as a general result. 

The linearized HPT-N3 theory is completed by inserting (|ll])-(|l3D into the linearized 
Euler equation (^. For the situations of interest here where all space dependence is in the 
z direction, this gives 

3 1 

= mnodtui - dzi-mvlnod^^ + -mvl^d^no) - hiFq - uqFi. (14) 



Here ^ is the fluid displacement in the z direction from 



i{z,t)= / ui{z,t')dt', (15) 



to 

and vp is a. position-dependent Fermi velocity determined by the inhomogeneous groundstate 
density riQ-. 

vl = ffm-\37r^no{z))^^\ (16) 

Equations ( p^lf ) - ([T^) are the basic equations of linear HPT-N3 hydrodynamicJ^Sl for a 
degenerate electron gas where both the groundstate and excitations exhibit one-dimensional 
spatial inhomogeneity. They were derived so as to satisfy the Harmonic Potential Theoremlli 
(HPT), to obey conservation of particle number and of momentum (hence Newton HI, N3), 
and to give the correct limit of high-frequency plasma dispersion in a uniform electron gas as 
in microscopic degenerate response (Lindhard-RPA) theory. We will show elsewhere^ that 
the HPT-N3 equations can alternatively be derived from nonlinear time dependent density 
functional theory with a very simple assumption regarding memory. The HPT-N3 equations 
represent an improvement over the HPT hydrodynamics originally proposedEl, in that the 
original form did not satisfy Newton HI. 

We will now outline the Tokatly and Pankratov (TP) hydrodynamic theoryil'lli, then 
compare it with the HPT-N3 theory summarized above, for the case of one-dimensional 
spatial variations and linear perturbations. Despite their very different derivations, these 
two theories will be shown to agree for this ID case. 

III. THE TP HYDRODYNAMIC MOMENT EQUATIONS FOR 
ONE-DIMENSIONAL SPATIAL VARIATION 

Tokatly and PankratovSS obtained hydrodynamic-style equations by taking velocity 
moments of a classical kinetic equation, assuming that this equation is adequate for degen- 
erate Fermi particles, presumably with a fully degenerate Fermi distribution in mind as the 
static solution. 
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They decoupled the set of hydrodynamic moment equations, using the smallness of the 
parameter 

7 = <'l^> (17) 

where < p/m > is the average particle speed, L = q^^ measures the spatial scale of the 
hydrodynamic perturbation, and Uc is an appropriate collision frequency. We are interested 
in the high-frequency, coUisionless case for which u » Uc and collision integrals can be 
ignored as discussed by TP. To 0(7^), their first three ID moment equations (Eqs. (16), 
(17) and (18) of their first paper@) are then, in our notation, 

d d d 
mi— + u—)n(z,t) + mn—u = (18) 
ot oz oz 

■'■ 'q^'"^ + ^('^'") ~ ^ {continuity) 
c) d d 

mn(— + u—)u + —P-nF = (19) 
ot oz oz 

{Newton II) 

{evolution of stress — pressure). 

Here the number density n, fluid velocity u, the stress-pressure tensor component P = P^z, 
and the total force F due to interactions, are all functions of the single space coordinate z 
and the time t. (The three-dimensional generalization of (0) - ( pO] ) is given by Equations 



(22) - (25) of TF'M, with L^^) and I^^) set to zero.) Our notation in (|l|)-(0) is chosen to 
match that of our recent wor The connections with the notation of TPlE3 follows, 
with our notation on the left: 



n = m 



-1^(0) 



V = -e0, 
F = +ed,<p. 

IV. COMPARISON OF LINEARIZED HPT-N3 AND TP THEORIES 

In linearizing ([ISD - (^Of) , Tokatly and Pankratov immediately replaced the convective 
derivative Dt = dt + u.V by the partial derivative dt- (See Eqs. (23)-(26) of TPlil and 
Eqs.(46) - (49) ff. of TP2§). In particular, this implies that quantities such as OPq/Oz and 
duo/dz SiTB set to zero, where Po{z) is the equilibrium pressure-stress tensor component and 
uqIz) is the groundstate density profile. Thus they specialized to systems with a spatially 
uniform groundstate, without actually stating this condition. They thereby also removed 
the possibility of treating edge inhomogeneities other than via assumed boundary condi- 
tions. Because of the homogeneous assumption they were furthermore unable to discuss the 
Harmonic Potential Theorem fully, even though some reference was made to it in Section 
V of TP2H. Note that the essential feature of the HPT, which makes it a more stringent 
test of many-particle theories compared with the Generalized Kohn Theorem^ (GKT), is 
that the entire groundstate many-body wavefunction, and hence all derived quantities such 
as no(r) and Po{r), are rigidly translated in HPT motion, giving n{r,t) = no(r — X{t)), 
P{r,t) = Po{r — X{t)) etc., where the center-of-mass coordinate X{t) satisfies the classical 
simple-harmonic-oscillator equation of motion. To verify satisfaction of the HPT, it is not 
sufficient to show that the centre of mass moves appropriately as in TP2lli: the preservation 
of the inhomogeneous groundstate spatial profiles no{r), Po{r) in the moving situation is also 
an essential feature. 
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A. Equilibrium state 



Fortunately it is easy to linearize (p!8[)-(pOD without discarding the space derivatives of 
groundstate quantities. First note that, for a 3D electron gas in a ID situation where 
spatial variation of all quantities occurs only in the z direction, the kinetic equation for the 
equilibrium distribution function /q is 

{-mr^VzQz + Fo.^}/o(^,p) = (21) 

where Fq = F^^z is the total selfconsistent force due to interactions. If -^0(2;) = 
— <9zVo(-2), then direct substitution into (|2l|) verifies that the equihbrium solutions are of 
the form 

where A{e) is an arbitrary function of one variable. For a degenerate electron gas, the 
classical kinetic theory will optimally decribe the quantal situation if A is chosen to be the 
Fermi distribution 

A{e) = Ce{ii-e), (22) 

where C and /i are independent of z and p. The groundstate density is then 

r rPFiz) 4^ 

no{z) = J fo{z,p)£p = C J Aixp^dp = C—pUz) 

where C is a constant and 

Pf{z) = mvp{z) = ^/2m{fi — Vo{z). 

Note that this assumption of a local Fermi distribution with a global chemical potential 
amounts to a Thomas-Fermi type of theory, which is therefore the natural groundstate 
theory to accompany a hydrodynamic treatment of excitations in a degenerate system. 

The zz component of the stress-pressure tensor from this distribution (see Eqs. (3) and 
(13) of TPlii) 
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is 



Pqzz = Lq^z = — [ PzPzfoir,p)(fp = I— I ffo{z,p)d-^p 
m 6m 



11/" 
Am J 



(23) 



Then 



z-^Qzz 



^z^Ozz 



dP( 



Ozz 



drio 



dzfiQ^^^ = d;,no^p%{z) = d^nj^vp{z) (24) 
6 no om 6 



B. Linearized moment equations of TP 



Linearizing our Eqs. (|T8D-(^) about the Thomas- Fermi groundstate just described, we 
write 

n(f, t) = no{z) + ni{z, t) 
u{f, t) = + ui{z, t) 

P,,{f,t) = Pozziz) + Piiz,t) 

F{f,t)=Foiz) + F,{z,t) 



and obtain 



Otn-i + no—- = 
oz 



mno—^ + d^Pi - noFi - niFo = 
ot 



(25) 
(26) 



dtPi + uid^Pozz + ^PozzdzUi = 



(27) 



where rii, ui and Pi are functions of z and t whereas no, -Foand Pozz are functions of z 
alone. Note that the uidzPozz term in (P7|) is missing in (e.g.) Eq. (49) of TP2 or Eq. 
(26) of TPl because TP have assumed a uniform unperturbed system so that Dt = df. It is 
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this additional term in our treatment which ensures satisfaction of the true HPT (not just 
the GKT which is satisfied by the treatment of TP involving uniform gas with zero-stress 
boundary condition). 

Insertion of the equilibrium pressure tensor and its space derivative from ( |23D and ( |2^ 
into ( P7| ) gives 

dtPi + ui{d,no)jvUz) + —vUz)nod,Ui = (28) 
-dtPi = —vl{z)uid^no + —vl{z)nod^Ui. (29) 

6 5 

We integrate both sides with respect to time, up to time t from an initial equilibrium at 
time to when all perturbations vanished. Then using (|15D we obtain 

Pi = -jvUz^no - ^vUz)nod,^. (30) 



Inserting (^0]) into (|26D we obtain 



dUi^ TTl 3777 

mno-^ - dz{—v],{z)id^no + —v],{z)nodzi) - uqFi - hiFq = (31) 

Eq. (^) derived from the Tokatly-Pankratov approach^ is identical to Eq. (|T^) derived 
from the HPT-N3 theoryS. This equivalence is the principal result of the present paper. 
Eq. (|3ll) has already been shownil to give a sensible treatment of both the Kohn mode and 
the hydrodynamic standing plasmon modes of the electron gas in a parabolic quantum well, 
while taking into account the smooth decay of the electron density at the edges and using 
only "natural" boundary conditions that the fluid displacement ^ and density perturbation 
rii are nowhere divergent. 

V. DISCUSSION 

We have compared two hydrodynamic theories describing plasmon excitations in degen- 
erate electron gases with one-dimensional spatial inhomogeneity. Their derivations differ 
principally by the way the pressure term is obtained. 
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The first theory, our HPT-N3 hydrodynamics, obtains the pressure term by identifying 
displacive and compressive components, and ni^, of the density perturbation rii. see 
Eqs. (0) - ([T0|). riiA is associated with the high-frequency pressure coefficient [3'^{uj oo) = 
^Vp, while is associated with the low-pressure coefficient = 0) = ^Vp. See Eq. 

(pup. (We will show elsewhere@ that these associations arise naturally from a form of time- 
delayed scalar local density approximation for the pressure.) By these means we were able 
to satisfy the Harmonic Potential Theoremlli (HPT) and the usual conservation laws, and to 
obtain within the same formalism the correct plasmon dispersion of the uniform electron gas. 
The HPT applies to harmonically confined systems and constitutes a stringent test of an 
inhomogeneous many-particle theory, requiring in the present case that there exist a Kohn 
mode in which the inhomogeneous groundstate density profile nQ{f) and pressure profile 
Po{r) move rigidly. This linear HPTN3 theory has been shown to give a sensible description^ 
not only of the Kohn mode and other (standing) modes of harmonically confined systems, 
but also of multipole surface plasmons on a near-neutral jellium electron gas layeiii. 

The second theory in our comparison is that of Tokatly and Pankratovil (TP) . This was 
derived for degenerate systems by assuming the validity of a classical kinetic equation, then 
truncating the momentum moment equations. This yields a prediction for the pressure term. 
The truncation is justified, in the case of interest to us, by the smallness of the parameter 
7 = vp/{Luj) where L and uj are the spatial scale and frequency. In order to compare with 
the HPT theory we needed to deal with inhomogeneous groundstates, and therefore we 
linearized the TP theory about the appropriate inhomogeneous Thomas-Fermi groundstate, 
a procedure not explicitly carried out by Tokatly and Pankratov. 

Despite their quite different derivations, these two linearized theories give the same 
predictions for the case of one-dimensional spatial variation, as evidenced by the identity of 
Eqs. (|T^ and (|3T|). This fact is the main result of the present paper. 

The HPT-N3 derivation shows that the assumption of a classical kinetic equation, made 
by TP even in degenerate cases, is not necessary. On the other hand, the TP derivation is 
part of a systematic expansion whereas the HPT-N3 derivation is not obviously part of any 
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systematic scheme. Thus the two approaches are somewhat complementary, and each tends 
to support the vahdity of the other. 

In order to be truly useful, the inhomogeous formalisms discussed here should ideally 
be applicable at edges where the groundstate and/or excited electron density may vary 
rapidly in space. It is then likely that the TP truncation parameter 7 (Eq. (|TTp) is small 
only at frequencies higher than the ones of interest. Nevertheless, in at least one important 
case covered by the Harmonic Potential Theorem, the formalism gives correct answers for 
very rapidly varying edge profiles. The use of these hydrodynamic approaches in regimes 
of rapid spatial variation is somewhat reminiscent of the commonplace and surprisingly 
successful use of the Local Density Functional formalismil for groundstate properties of 
highly inhomogeneous systems, despite the in-principle restriction to slow spatial variation. 
In that case the success is at least partly explained by the satisfaction of sum rules and 
constraints^. The HPT constraint can be viewed in the same light for time-dependent 
cases, and indeed has already been used for this purpose in the context of the time- dependent 
Local Density Approximation for exchange and correlation^. The connection of the present 
hydrodynamic approximations with finite-memory versions of density functional theory will 
be made more explicit elsewhere^. 

It should also be stressed that the existing HPT-N3 theory does not necessarily agree with 
the more general three-dimensional versionEl of the TP theory. This is because the HPT-N3 
theory to date has assumed a scalar pressure, whereas in the high-frequency plasmon case 
the pressure should certainly be a tensor, which indeed is what emerges from the 3D TP 
theory. It will be interesting to see if a tensor ansatz for the pressure, along the lines of the 
HPT-N3 argument, can re-derive the 3D TP theory. 
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